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ABSTRACT 

We analyze the three catalogs of nearby loose groups by Garcia (1993). She identified groups in a 
magnitude-limited redshift galaxy catalog, which covers about ~ 2/3 of sky within cz = 5500 km s~^, by 
using two methods, a percolation and a hierarchical method. The free parameters of the group-selection 
algorithms were tuned to obtain similar catalogs of groups. The author also proposed a third catalog of 
groups defined as a combination of the two. Each catalog contains almost 500 groups. 

In agreement with previous works on earlier catalogs, we find that groups can be described as collapsing 
systems. Their sampled size is in general considerably larger than their expected virialized region. We 
compute the virial masses and correct them by taking into account the young dynamical status of these 
groups. We estimate corrected group masses, M, for two reference cosmological models, a flat one with 
a matter density parameter f2o = 1 and an open one with Qq = 0.2. For each of the three catalogs we 
calculate the mass function. 

We find that the amplitude of the mass function is not very sensitive to the choice of the group- 
identification algorithm. The number density of groups with M > 9 x 10"'^^ M©, which is the 
adopted limit of sample completeness, ranges in the interval 1.3-1.9 x lO^^h^Mpc'^ for ^lo = 1, and it is 
about a factor of 15% lower for = 0.2. The mass functions of the hierarchical and combined catalogs 
have essentially the same shape, while the mass function of the percolation catalog shows a flattening 
towards large masses. However, the difference decreases if we do not consider the most massive groups, 
for which reliable results come from galaxy cluster studies. 

After having estimated the mass contained within the central, presumably virialized, regions of groups 
by adopting a reduction in mass of ~ 30-40%, we do a comparison with the results coming from the 
virial analysis of nearby rich clusters (Girardi et al. 1998a). All three group mass functions turn out 
to be a smooth extrapolation of the cluster mass function at M < 4 x 10^'^h~^ M©, which is the 
completeness limit of the cluster sample. The resulting optical virial mass function of galaxy systems, 
which extends over two orders of magnitude, is fitted to a Schechter expression with a slope of ^ —1.5 
and a characteristic mass of M* ~ 3 x 10^'^h~^ Mq. We also verify that our group mass function 
reasonably agrees with the Press-Schechter predictions of models which at large masses describe the 
virial mass function of clusters. 

Subject headings: galaxies: clusters: general - cosmology: observations - cosmology: theory - large scale 
structure of universe. 
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small number of group members and to uncertainties in the 
dynamical stage. In fact, although group cores are virial- 
ized or close to virialization (Zabludoff & Mulchaey 1998), 
the size of groups identified in three-dimensional galaxy 
catalogs, i.e. ~ 0.5-1 Mpc, is appreciably greater than 
their expected virialized region, i.e. ^ 0.2—0.4 Mpc for 
systems with a line-of-sight velocity dispersion of 100-200 
km (according to the relations found for galaxy clus- 
ters, e.g. Carlberg, Yee, & Ellingson 1997; Girardi et al. 
1998b). Indeed, there is a strong indication that these 
groups are not virialized systems over the whole sampled 
region, but can be rather described as being in a phase of 
collapse (e.g., Giuricin et al. 1988; Mamon 1994; Diaferio 
et al. 1993). Therefore, usual estimates of velocity disper- 
sion and virial mass are not easily connected to physical 
quantities such as group potential and mass. 

The small number of data and the uncertainties on dy- 
namical status prevent one to use refined methods to re- 
ject interlopers in each individual group (e.g.. Zabludoff & 
Mulchaey 1998, Mahdavi et al. 1999). Instead, one must 
rely on member galaxies as assigned by the group-selection 
algorithm, while checking a posteriori the presence of spu- 
rious groups in a statistical sense (e.g., Ramella, Pisani, 
& GoUcr 1997; Diaferio ot al. 1999). Indeed the results 
could depend on the choice of the group-sekx'tion algo- 
rithm and its free parameter (e.g., Pisani et al. 1992 - 
hereafter P92; RamcUa ct al. 1997). For instance, Fred- 
eric' s (1995b) analysis of cosmological N-body simulations 
suggested that the estimated median mass depends on the 
algorithm and that the resulting bias is sensitive to the 
depth of the galaxy survey. However, even the analysis 
of simulated groups is not an easy task and, indeed, the 
results on mass can depend on the treatment of halos (cf. 
Frederic 1995b). 

A further uncertainty is connected to cosmic variance. 
In fact, group catalogs are recovered from local galaxy cat- 
alogs which may not be fair samples of the universe. 

In view of these difficulties, few statistical distributions 
of group dynamical properties are available in the litera- 
ture and they are often discrepant. The cumulative dis- 
tributions of internal velocity-dispersion, as computed by 
Moore, Frcnk, & White (1993) and by Zabludoff ct al. 
(1993), are strongly discrepant (the number densities of 
groups with line-of-sight velocity dispersion larger than 
200 km s^^ differ by a factor of 100, see Fig. 6 of Fadda 
et al. 1996 for a comparison). Moreover, analyzing nearby 
groups {cz < 2000 km s~^) of three different group cata- 
logs, P92 found a significant dependence of the distribu- 
tion of mass and other dynamical group parameters on the 
group-identification algorithm. 

The availability of new group catalogs has prompted 
us to derive a new group mass function, whose connec- 
tion with the recent determination of the optical virial 
mass function of nearby rich galaxy clusters (Girardi et al. 
1998a, hereafter G98) deserves to be investigated. 

The work by Garcia (1993, hereafter G93), who con- 
structed two group catalogs using two difFc;rcnt identifica- 
tion algorithms (the percolation and hierarchical ones) and 



proposed a third catalog which is a combination of the two, 
represents a good data base for facing the effect of iden- 
tification algorithms. So far, G93 catalogs are the largest 
catalogs of groups presently published. They are largely 
superior to those analyzed by P92 both for the number 
of groups (450-500 groups for each of the three catalogs) 
and the encompassed volumes (^ 2/3 of sky, cz < 5500 
km s~^). Moreover, these group catalogs were selected 
from the same parent galaxy sample, thus allowing us to 
better investigate on the effects due to differences in the se- 
lection algorithm. Furthermore, the improved statistics in 
the high-mass range (less than ten groups analyzed by P92 
have masses larger than 10^'^h~^ Mq) permits an interest- 
ing comparison with cluster mass function and a determi- 
nation of the virial mass function over an unprecedently 
large range of masses. 

In § 2 we briefly describe the data. In § 3 we calcu- 
late group masses. In § 4 and 5 we compute group mass 
function and verify its stability, respectively. In § 6 we 
compare the results of groups and clusters, recovering the 
mass function of galaxy systems for a mass range which 
extends over two orders of magnitude. In § 7 wc give our 
discussions. In § 8 we summarize our results and draw our 
conclusions. 

Throughout the paper, errors are given at the 68% confi- 
dence level and the Hubble constant is Hq = 100 h Mpc"-^ 
km s~^. 

2 DATA SAMPLE 

We analyze the loose group catalogs constructed by 
G93. These groups were identified by using galaxies 
(within cz = 5500 km s~^) belonging to the subsam- 
ple of the Lyon-Meudon Extragalactic Database which is 
nearly complete down to the limiting apparent magnitude 
Bq = 14, the total blue magnitude corrected for Galac- 
tic absorption, internal absorption and K-dimming. G93 
used two methods in group construction: a percolation 
method (hereafter P, derived from the friends of friends 
method presented by Huchra & Geller 1982) and a hier- 
archical method (hereafter H, derived from that of Tully 
1980, 1987). Each method gives one catalog. The P and 
H catalogs contain 453 and 498 groups of at least three 
members, respectively. In particular, G93 tuned the free 
parameters of the methods so as to obtain the best com- 
promise between the stability of group membership and 
the similarity of the two group catalogs (Garcia, Morenas, 
& Paturel et al. 1992). 

Then G93 combined together the two catalogs in or- 
der to obtain the final catalog (hereafter G catalog, 485 
groups) defined as the catalog which contains only the 
groups which were found in part in both catalogs. For 
these final groups only the galaxies in common were kept 
as group members. If some groups of a catalog (in most 
cases the P catalog) turned out to be divided into two or 
more groups in the other catalog (in most cases the H cat- 
alog), the smallest systems were kept for the final catalog. 

In our work, as already suggested by the author, we do 
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not consider galaxies added afterwards (flagged by "+") 
which have no known magnitude or which are not fuhy 
satisfying the selection criteria. We refer to the original 
paper for more details on group catalogs. 

We exclude from our analysis groups with cz < 500 
km because, when the velocity becomes low, its ran- 
dom component dominates and the velocity is no longer a 
reliable indication of the distance. 

The samples we analyze consist of 446, 490, and 476 
groups (P, H, and G groups, respectively). 

3 COMPUTING GROUP MASSES 

The virial mass of a relaxed galaxy system is computed 
as — (T^Rv/G where Ry is the virial radius of the 

system and a is the velocity dispersion of member galax- 
ies (Limber & Mathews 1960). Other usual mass estima- 
tors, e.g. the projected and the median mass mass esti- 
mators (Bahcall & Tremaine 1981; Heisler, Tremaine, & 
Bahcall 1985), give similar results as shown by P92 (cf. 
also Pcrca, del Olmo, & Moles 1990). A more serious 
problem comes from the fact that one must assume that, 
within each group, mass distribution follows galaxy distri- 
bution (e.g., Mcrritt 1987). This assumption is shown to 
be enough reliable for galaxy clusters both from optical. 
X-ray, and gravitational lensing analyses (see Girardi et 
al. 1998b, Lewis ct al. 1999 and references therein) and, 
for likeness, the same assumption can be made for galaxy 
groups. 

We compute the above virial parameters from the ob- 
served projected positions in the sky and the line-of- 
sight velocities of the member galaxies. In fact, for a 
spherical system, the parameters Ry and a are linked 
to their observational counterparts as a = V^a^ and 
Rv = {ir/2)Rpv = {TT/2)Nm{Nm - l)/^i>jR-j\ where 
(T„ is the line-of-sight velocity dispersion, Rpv is the pro- 
jected virial radius, is the number of group mem- 
bers and Rij the galaxy projected distances. In particu- 
lar, we estimate the "robust" velocity dispersion by using 
the biweight estimator for rich groups (member number 
Nm > 15) and the gapper estimator for poorer groups 
(ROSTAT routines by Beers, Flynn, & Gebhardt 1990). 
Beers et al. (1990) have shown the superiority of their 
techniques in terms of efficiency and stability when treat- 
ing systems with a small number of members (cf. also 
Girardi et al. 1993). In our case (see also Mahdavi et 
al. 1999) we verify that the distributions of robust and 
traditional estimates of velocity dispersion are not dif- 
ferent according to the Kolmogorov-Smirnov test (here- 

TABLE 1 

Group Properties 



after KS-test, e.g. Ledermann 1982). We apply the rel- 
ativistic correction and the usual correction for velocity 
errors (Danese, De Zotti, & di Tullio 1980). In partic- 
ular, for each galaxy, we assume a typical velocity error 
of 30 km based on the average of errors estimated 
in RC3 catalog from optical and radio spectroscopy (de 
Vaucouleurs et al. 1991). 

Note that our estimates of the virial parameters do 
not require any luminosity-weighting procedure. In- 
deed, it was shown that the values of the virial masses 
are largely insensitive to different weighting procedures 
(e.g., Giuricin, Mardirossian, & Mezzetti 1982; P92) in 
agreement with evidences of velocity equipartition (e.g., 
Giuricin et al. 1982). 

In order to take into account the dynamical state of 
groups we use the method proposed by Giuricin et al. 
(1988). This method is based on the classical model of 
the spherical collapse where the initial density fluctuation 
grows, lagging behind the cosmic expansion when it breaks 
away from the Hubble flow, and begins to collapse; then 
a relaxation process sets in (e.g., Gott & Rees 1975). The 
time evolution curves A{t), which is the; ratio between the 
absolute values of the kinetic and potential energies, given 
by the dynamical model, are the starting point for the 
determination of the evolutionary stage. The evolution 
curve used by Giuricin et al. (1988) has been derived from 
numerical simulations of systems composed of 15 point 
masses with a Schechter-type mass function (Giuricin et 
al. 1984; the limits of this model are discussed in §7.1) 

According to the above method, the value of A, 
which is needed to recover corrected masses as M = 
(l/2A)((T^i?y/G), can be inferred from the estimate 
of the presently observed virial c;rossing time tcr = 
(3/5)"^/^i?y/cr. One can also derive the value of t, which is 
the elapsed time since fluctuations started growing (here 
in units of the crossing time at the virialization, t"^). In 
particular, the precise values of A and r depend on the 
background cosmology. Here, we estimate corrected group 
masses for two reference cosmological models, a flat one 
with = 1 for the matter density parameter and an 
open one with fig = 0.2 (cf. also P92). 

For each catalog, in Table 1 we give the number of 
groups (Col. 2) and median group properties: the mean 
redshift (Col. 3); the line-of-sight velocity dispersion, cfy 
(Col. 4); the projected virial radius, Rpv (Col. 5); the 
crossing time, tcr (Col. 6); the virial mass, A^vir (Col. 7); 
the value of r and the corrected virial mass, M, for the 
ilo = 0.2 model (Cols. 8 and 9) and for the fio = 1 model 
(Cols. 10 and 11). 
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From the values of t (< 2tt) we infer that most groups 
are in the phase of coUapse and not yet viriahzed in 
agreement with previous analyses of earlier catalogs (e.g., 
Giuricin et al. 1988; P92). Moreover, P groups result to 
be more evolved than H groups. The resulting estimate 
of the corrected mass is larger than the virial mass by a 
factor of 20-60%, depending on the group catalog and on 
the assumed background cosmology. Since, for each of the 
three catalogs, tcr does not correlate with Myir, one could 
apply the same percent correction to M-^ir of all groups 
to obtain corrected masses (as inferred from Table 1, e.g. 
multiplying by 12.7/8.4 = 1.41 in the case of G groups 
and rio = 1). The mass distribution resulting from the 
application of an average correction is indistinguishable 
(according to the KS-test) from the application of individ- 
ual corrections for each group. 

The distribution of corrected group masses is only 
slightly dependent on the cosmological environment; in 
fact, the two mass distributions computed for fio = 0.2 
and fio = 1 do not differ significantly (according to the 
KS-test). Therefore, hereafter, if not explicitly said, we 
consider the f2o = 1 case only. 

4 GROUP MASS FUNCTION 

In order to compute a reliable group mass function, MF, 
we avoid strongly obscured regions by considering only 
groups with galactic latitude |6| > 20° (i.e. within a solid 
angle of a; = 8.27 sr). We consider 409, 381, and 344 G, 
H, and P groups, respectively. 

To obtain the true spatial number density, the selection 
effects of the galaxy catalog, which is magnitude limited, 
must be taken into account. Following Moore et al. (1993), 
we weight each group by using the magnitude of its third 
brightest galaxy which allows the inclusion of the group it- 
self in the catalog. We weight each group by w = 1/Tmax 
where Tmax is the maximum volume, within the catalog 
velocity cut-off, out to which the group can be seen: 
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where vum is the smaller one between 5500 km and the 
maximum recession velocity at which the third brightest 
galaxy in the group would be brighter than the magni- 
tude limit; c is the speed of light and g'o the deceleration 
parameter. 

By plotting the weights vs. group masses (cf. Figure 1, 
left panels), we note that groups generally lie in a well- 
defined region of the w — M plane with the exception of a 
few groups which fall very far away from the other points. 
Therefore, in order to avoid problems of instability in the 
resulting MF, for a few groups we recompute the weights 
according to the following procedure. In each mass range 
(of a unity in logarithmic scale) we compute the mean and 
s.d. of logw and substitute the values which are three s.d. 



far away from the mean with the corresponding mean val- 
ues. We change the value of weights for three, four, and 
three groups in G, H, and P catalogs, respectively (cf. 
Figure 1, right panels). 
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Fig. 1. — For the three catalogs, left panels show that groups 
generally lie in a well-defined region of the weight-mass plane. 
Right panels show the same, after recomputing the value of 
weights for a few groups which lie far from this region (see 
text). 

Then, since the parent galaxy sample was foTind to have 
aredshift incompleteness of ~ 10% (Marinoni et al. 1999), 
the resulting group densities are enhanced by the same fac- 
tor. This is a rough correction. As discussed in § 5, a more 
refined correction which takes into account the number of 
group members would give similar results. 

The comparison between different catalogs is shown in 
Figure 2. In the low-mass range the three MFs show 
large fluctuations (well over the estimated Poissonian er- 
rors) and a tendency towards flattening, which suggests 
problems of incompleteness. Our results for M ^ 9 x 10^^ 
Mq, which we assume as our completeness limit, come 
from 230, 207, and 214 groups for G, H, and P catalogs, re- 
spectively. The determination of global group number den- 
sity is quite stable, since the density of groups with M > 
9 X 10i2/;-i A4q is 1.4, 1.3, and 1.9 x lO'^ {h''^ Mpc)'^ for 
G, H, and P catalogs, respectively. 

However, Figure 2 shows that P catalog gives a flatter 
MF with respect to G and H catalogs. In fact, according 
to the KS-test, the MF of P groups strongly differs from 
that of G and H groups (at a c.l. greater than 98%), while 
there is no difference between G and H groups. However, 
the difference becomes smaller if we exclude large masses 
and, for instance, we do not find any significant diff'erence 
if we exclude M > 4 x lO^^ft-^^ Mq, a range where reliable 
results come from cluster analysis. 
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In the case of the fta = 0.2 model wc find similar re- 
sults, except for the fact that, since estimated masses are 
smaller, the density of groups with M > 9 x 10^^h~^ Mq is 
somewhat lower than in the l^o = 1 case, i.e. 1.2, 1.0, and 
1.7 X 10~^ {h~^Mpc)~^ for G, H, and P catalogs, respec- 
tively. 
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Fig. 2. — We compare mass functions as computed from the 
three catalogs. G, H, and P groups are denoted by closed 
circles, open circles, and triangles, respectively. To avoid con- 
fusion we give error bars {la Poissonian uncertainties) only for 
G and P groups. 



5 ON THE STABILITY 
OF GROUP MASS FUNCTION 
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Fig. 3. — The cumulative mass function for the nearby groups 
of G, H, and P catalogs (with cz < 2000 km ). Error bars 
represent la Poissonian uncertainties. 

First we consider groups with cz < 2000 km s~^. In 
particular, we have 136, 124, and 117 G, H, and P nearby 

groups, respectively, with galactic latitude |6| > 20°. For 
each catalog, the observed group population in the nearby 



subsamples can be considered a good representation of 
the total population because we do not see any signif- 
icant trend between mass and distance. Therefore, the 
nearby groups can be considered a complete sample, ex- 
cept for the less massive ones, which could not be iden- 
tified out to the distance limit of 20 Mpc since their 
third brightest galaxy is not bright enough. If we assume 
that the nearby groups form a complete sample, we can 
directly compute the MF for each of the three catalogs. 
There is no significant difference among the three mass 
distributions (according to the Kruskas Wallis test, e.g., 
Ledermann 1982, cf. Figure 3). Moreover, apart from the 
density of low-mass groups, M ^ 10^^ h^^ Mq, where the 
completeness of the nearby groups is questionable, there is 
a good agreement between the number densities estimated 
from the whole catalogs with those from the nearby ones 
(cf. Figure 4). This result reassures us of the reliability of 
our weighting procedure (cf. § 4). 



■* 
o 



a 
a 



0.1 
0.01 
0.001 
0.0001 

10-5 
10-6 
10-^ 

0.1 
0.01 
0.001 
0.0001 

10-5 

10-6 
10-^ 

0.1 
0.01 
0.001 
0.0001 
10-5 

10-6 

10-'^ 



G 



I I llllllj I I llllllj 



I I lllllll I I II 



via 

■ n 



+H 



■ ^ — 



A a 
! n 



r H 



■Hi 



^ (D* • 



ays 



P 



iiiiiil 



I III 



jL 



I I lllllll 



ill 



1012 1013 1014 1015 

M / h-1 

Fig. 4. — For the three catalogs we compare mass functions as 

computed from all the grotips and the nearby groups (closed 
and open circles, respectively). Error bars represent la Pois- 
sonian uncertainties. 

The physical reality of the detected groups is often dis- 
cussed in the literature. In particular, the efficiency of 
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the percolation algorithm lias been repeatedly checked 
through cosmological N-body simulations (e.g., Nolthe- 
nius & White 1987; Moore et al. 1993; Nolthenius, Klypin, 
& Primack 1994; Frederic 1995a, b; Nolthenius, Klypin, 
& Primack 1997; Diaferio et al. 1999) and geometrical 
Monte-Carlo simulations (Ramella et al. 1997). These 
computations show that an appreciable fraction of the 
poorer groups, those with Nm < 5 members, is false (i.e. 
unbound density fluctuations), whereas the richer groups 
almost always correspond to real systems (e.g., Ramella et 
al. 1995; Mahdavi et al. 1997). By following the results of 
Ramella et al. (1997; cf. also Diaferio et al. 1999) we re- 
duce the weights of groups with Nm = 3 and 4 < Nm < 5 
members by 70% and 20%, respectively. Figure 5 shows 
the effect on the cumulative MF of P groups. The MF 
at large masses is quite stable and, however, the number 
density at the completeness limit of 9 x 10^^ h^^ Mq could 
be overestimated only by a factor of ^ 50%. 




1012 IQ13 IQ15 
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Fig. 5. — We show the effect of the presence of poor groups 
on the cumulative mass function of P groups. The thin line 
indicates the mass function obtained when the weight of poor 
groups is suitably reduced (see text). Error bars represent la 
Poissonian uncertainties. The mass completeness limit is in- 
dicated by the vertical dashed line (M > 9 x 10^^ Mq). 

Galaxy density is known to show significant fluctua- 
tions around the mean density. Redshift surveys reveal 
voids of sizes up to 50 Mpc and large bidimensional 
sheets (e.g., de Lapparent, Geller. & Huchra 1986; Geller 
& Huchra et al. 1989; Vettolani et al. 1997). Indeed, 
there is also an indication of a local underdensity of the 
above size (e.g., Zucca et al. 1997; Marinoni et al. 1999). 
Since group catalogs are obtained from galaxy catalogs, 
one could suspect that the amplitude of the group MF we 
estimate in the local universe is far from being a fair value. 
In order to shed light on this point, we use the groups iden- 
tified in the northern Center for Astrophysics redshift sur- 
vey (CfA2N) by Ramella et al. (1997). The CfA2N covers 
a smaller solid angle (1.2 sr) than the survey analyzed by 
G93, but it is deeper {cz < 12000 km s~^) and contains the 
Great Wall, a very overdense structure . From Ramella et 
al. (1997) we take group dynamical quantities, Myir and 



tcr, and apply the same procedure outlined in § 3 and 4 
for the evolutionary correction and for the computation 
of weights (we take the luminosity of the third brightest 
galaxies from Ramella et al. 2000, in preparation). The 
resulting MF is shown in Figure 6. The number density of 
the CfA2N groups with M > 9 x lO^^/j-i ij^g within 
the range of values we found for G93 catalogs. In con- 
clusion, although the CfA2N groups come from a quite 
different volume, their MF is similar to the MF of the 
G93 groups. This result agrees with that by Ramella et 
al. (1999) who found that different group catalogs, which 
sample different volumes, give similar velocity-dispersion 
distributions. 
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Fig. 6. — We compare the cumulative mass function of the 

CfA2N groups by Ramella et al. (1997; RPG97) with that 
of G, H, and P groups. Error bars represent la Poissonian 
uncertainties. The mass completeness limit is indicated by the 
vertical dashed line (M > 9 x 10^^ Mq). 

In § 4 in order to take into account the 10% redshift 

incompleteness of the parent galaxy sample (Marinoni et 
al. 1999), we have applied a rough, very small upward cor- 
rection by 10% to the group densities. A more refined cor- 
rection should take into account the number of members 
in each group. For instance, let's adopt the extreme view 
that the incompleteness of parent galaxy sample does not 
affect the density of groups with more than three members, 
but affects only groups with three members. In the worst 
case, out of 100 groups (i.e. 300 galaxies, of which 30 are 
missed) one misses 30 groups. According to this kind of in- 
completeness we recompute the MFs for the three catalogs. 
The differences in densities for groups with M > 9 x 10^^ 
are very negligible for G groups and at most amount to 
<^ 10%, which is within the errors, for H and P groups. 

6 GROUP VS. CLUSTER MASS FUNCTION 

Contrary to the case of groups, recent determinations of 
the distribution of internal velocity dispersions for nearby 
rich clusters well agree within the errors (cf. Figures 6 
and 9 of Fadda et al. 1996). The cluster MF based on 
optical virial masses has recently been presented by G98. 
In particular, masses were computed within a radius en- 
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closing the region where chisters are virialized. The re- 
sulting MF is reliably estimated for masses larger than 
4 X 10^^h~^ Mq with ~ 50 nearby clusters available for 
this mass range. 

We must take into account that the groups examined 
here extend outside their virialized regions for a meaning- 
ful comparison of group masses with cluster masses as well 
as with Press & Schechter (1974, hereafter PS) predictions, 
which hold for virialized objects (e.g., Eke, Cole, & Prenk 
1996; Laccy & Cole 1996; Borgani et al. 1999). 

By using eq. 3 of Ramella et al. (1997), we estimate 
that P groups has a limiting number density contrast of 
~ 70, while for H groups, detected imposing a luminosity 
density threshold (Gourgoulhon, Chamaraux, & Fouque 
1992; G93), we compute a luminosity density contrast of 
~ 40. These calculations use the Schechter luminosity 
function (1976) with a slope of —1.1, a normalization fac- 
tor of 0.014 Mpc^^, and a characteristic magnitude of 
Mg = —20.0 — 5 • logh, as well as the galaxy blue lumi- 
nosity density, ^ 4.5 x 10® Lq h'^ Mpc^^, as obtained by 
Marinoni et al. (1999) for the same parent galaxy sample. 
The estimated values of density contrast should be consid- 
ered as rough estimates also in view of the difficulties in 
obtaining group catalogs with a constant density contrast 
(e.g., Nolthenius & White 1987). 

In order to compute the matter density contrast, Sp/p, 
from the above galaxy density contrast, {6p/p)g, one 
should know the biasing factor b = {6p/p)g/{6p/p). Here- 
below we compute the value of bias as 6 = l/crg where 
(78, the r.m.s. mass density fluctuation in spheres of 8 
Mpc radius, is estimated for different values of flo ac- 
cording to the G98 relation found by comparing PS pre- 
dictions to the cluster MF (cf. their cq. 4). Accordingly, 
we adopt h = 1/0.60 ~ 1.7 and b = 1/1.23 ~ 0.8 in the 
case of the Qq = 1 and the = 0.2 models, respectively. 

The resulting Sp/p for P and H groups (41 and 24, re- 
spectively, for = 1; 88 and 50 for = 0.2) are much 
smaller than the values of ~ 180 and ^ 550 expected 
within the virialized region for fig = 1 and CIq = 0.2, 
respectively (e.g., Eke et al. 1996). 

After assuming that groups have a common radial profile 
(Fasano et al. 1993), we can roughly estimate the fraction 
of the mass contained in the virialized region. For each 
of the three catalogs. Figure 7 plots the cumulative distri- 
butions of the galaxy distances from the group (biweight) 
center combining together data of all groups. To combine 
the galaxies of all groups we divide each galaxy distance to 
the projected virial radius, Rpv, of its group. Moreover, 
in order to better show the behavior of galaxy density, we 
also normalize the distances to the mean < Rmax/Rpv > 
of the catalog, where the maximum radius. Rmax, is the 
projected distance from the group center of the most dis- 
tant galaxy, and we normalize the numbers of objects. A*", 
to that c;ontained within < R,nax/Rpv >■ 

From Figure 7 one can infer the fraction of the num- 
ber of galaxies (i.e. the fraction of group mass if galaxies 
trace mass) contained within each radius (i.e. how density 
changes with radius). We are interested in determining 



the radius (and the corresponding galaxy number/group 
mass) for which one obtains a density enhancement which 
is large enough to reach the density contrast expected at 
virialization. In particular, in the case of fig = 1 a reduc- 
tion in radius of ~ 50% (i.e. in number/mass of ~ 30%) is 
enough to reach the virialization density. In fact, from Fig- 
ure 7 one infers that 70% of galaxies are contained within 
about half of the radius (0.545 for P groups; 0.49 for H 
groups). Therefore, considering only these central group 
regions, the density increase is = 0.7/(0.545'^) ^ 4.3 for 
P groups and 0.7/(0.49*^) ~ 5.9 for H groups; in other 
words, the resulting density contrast in a = 1 universe 
is ~ 41 X 4.3 ~ 180 for P groups and 24 x 5.9 ~ 140 for 
G groups (these values are comparable to the virialization 
density contrast). In view of the inherent approximations, 
we simply adopt a reduction of group masses by 30% for 
all groups in each of the three catalogs. In the case of 
= 0.2, similar arguments lead us to adopt a reduction 
in radius of ^ 60% (i.e. in number of ^ 40%) to reach the 
virialization density. 

In the following comparison with clusters we will reduce 
group masses by 30% and 40% for the ilo = 1 and the 
f2o = 0.2 cases, respectively. We remark that this kind 
of correction is taken to be criual for all groups since we 
assume a common typical spatial distribution of member 
galaxies (Fasano et al. 1993) and we do not consider each 
individual evolutionary state (note that the crossing time 
does not show any correlation with group mass, cf. the 
end of § 3). 
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(R/Rpv)/<Rpv/R,„ax> 

Fig. 7. — For each of the three catalogs we give the cumulative 
distribution of the galaxy distances from group center. To com- 
bine the galajcies of all groups we divide each galaxy distance 
to the projected virial radius, Rpv, of its group. Moreover, we 
also normalize the distances to the mean < Rmax/Rpv > of 
the catalog, and normahzc the numbers of objects, N, to that 
contained within < Rmax/Rpv >■ The distributions of G and 
H are nearly overlapping. 

Figure 8 shows the group and the cluster MFs. The P 
group MF gives too high values at high masses when com- 
pared with cluster results, while G and H catalogs do not 
contain very massive groups. Therefore, the differences 
among group catalogs in describing the MF at high masses 
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(cf. § 4) lead to strong uncertainties just where cluster 
data are available. However, for M < 4 x 10^'^h~^ Mq all 
the three group MFs can be regarded as a smooth extrap- 
olation of the cluster MF. 

To provide a phenomenological fitting to the data which 
takes into account the effect of mass uncertainties we use 
the Schechter expression (1976). A theoretical, PS ap- 
proach is instead discussed at the end of this section. 

Following a maximum-likelihood approach we fit the 
Schechter expression for our MF on the whole mass range 
(groups with 9 x W^h'^ Mq < M < 4x lO" h'^ Mq and 
clusters with M > 4 x lO" h'^ Mq): 



3.9 X lO^^/i-^ Mp,. 



n{M) = n* 



M 



-M/M' 



(2) 



where n(M) is suitably convolved with the mass errors 
AM. We compute uncertainties on group masses in the 
same way as in Girardi et al. (1998b): we propagate statis- 
tical errors in the estimate of the velocity dispersion and 
of the virial radius Rpv , for which errors are estimated via 
bootstrap and jacknife techniques, respectively. The re- 
sulting mass uncertainties range from 15% for clusters at 
M ~ 2 X 10^5 /i-i M© to - 90% for groups at M ~ 9 x lO^^ 
h-^ Mq (we fit log{AM/M) = 4.548 - 0.355 • log M) and 
we assume a lognormal error distribution. The effect of 
these uncertainties is negligible in the cluster mass range, 
but it becomes significant in the range of low mass groups. 

For instance, in the case of G groups (Oq = 1 model) 
the resulting fitted parameters are: n* = (2.2 ± 0.3) x 
10-5 (/i-iMpc)-3 (1014/i-iMq)-i, M* = S.ltli X 10" 
/i"^ Mq, and a = 1.55to\t- Table 2 gives the fit results 
for the three catalogs and for both cosmological model: 
(Col. 3) gives the number of groups used in this analysis; 
(Cols. 4, 5, and 6) give the fitted n*, M*, and a param- 
eters, respectively. The values of a vary in the range of 
1.3-1.6 and the values of M* vary in the range of 2.7- 
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Fig. 8. — For the two reference cosmological models and for 
each of the three catalogs we show group and mass functions, 
where masses are computed within the virialized region. Clus- 
ter data come from Girardi et al. (1998a). We show the fitted 
Schechter function and its convolution with errors (thin and 
thick lines, respectively) by combining data of G groups and 
clusters. 



TABLE 2 

Parameters of the Mass Function Fit 



Cat. 


Model 


N 


n* 


M* 


a 








[{/i-iMpc)-3 (10" h-i Mq)-1] 


(lO"h-iM0) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


G 


no = 1 


193 


2.2x10-5 


3.1 


1.55 


H 


no = 1 


177 


1.0x10-5 


3.9 


1.64 


P 


no = 1 


153 


4.6x10-5 


3.1 


1.42 


G 


no = 0.2 


159 


2.4x10-5 


2.7 


1.49 


H 


no = 0.2 


147 


1.5x10-5 


3.1 


1.55 


G 


no = 0.2 


153 


5.5x10-5 


3.0 


1.31 



The errors on the "shape" parameters a and M* are 
directly given by the maximum likelihood method (Avni 
1976). We recover the error on the normalization n* by 
considering the Poissonian error bars associated to the 
global number of objects considered (and assuming the 
best fit parameters for the shape). These errors are the for- 
mal ones and do not consider other additional effects. The 
effect of incompleteness (see the end of § 5) is very small. 



being smaller than Poissonian errors. The effect due to the 
presence of interlopers is evidenced by differences in results 
relative to different group-selection algorithms. In particu- 
lar, for P groups, for which the presence of spurious groups 
is well studied in the literature, adopting the correction 
analyzed in § 5, we find values for the slope and charac- 
teristic mass (a = 1.5 and M* = 3.3 x lO^ft.-^ Mq for 
the case flo = 1, respectively) which are within the errors, 
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but a n* smaller by a factor of 50%, e.g. n* = 2.4 x 10~^ 
{h~^Mpc)^'^ (lO^'^/i^^M©)^^ as already pointed out in 
§ 5. However, this normalization is similar to that coming 
from the G catalog. 

G98 applied the PS approach to constrain the cosmolog- 
ical parameters from the observational cluster mass func- 
tion. The PS approach provides fairly accurate analytical 
approximation to the number density of dark matter halos 
of a given mass. The halo mass which appears in the PS 
formula refers to the mass contained within the virialized 
region, i.e. the region with a present density contrast of 
IStt^ = 178 (for Oq = 1, see Eke et al. 1996 for values 
in different cosmologies). Moreover, the PS mass function 
has been compared with N-body simulations by several 
authors (e.g.. Eke et al. 1996; Lacey & Cole 1996; Gross 
et al. 1998; Borgani ct al. 1999; Governato et al. 1999) 
and has been generally shown to provide a rather accurate 
description of the abundance of virialized halos of cluster 
size. 
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Fig. 9. — For the two reference cosmological models and for 
each of the three catalogs we compare observational mass func- 
tions with Press-Schechter mass functions, with and without 
convolving with the uncertainties in the mass estimates (thick 
and thin lines, respectively) . The plotted Press-Schechtcr mass 
functions are those found to well describe cluster data (Girardi 
et al. 1998a). 

It is worth verifying that the present group MF, in spite 
of all difficulties inherent in the analysis of groups (cf. 



§7.1), is a good extension of the PS form fitted on clusters. 
We use the same PS approach used by G98 for describing 
the MF of galaxy clusters (cf. their eq. 2 and 3 and details 
in § 3). In particular, G98 recovered the relation between 
as and (cf. their eq. 4 for f^A = 0). Accordingly, we as- 
sume as = 0.60 and = 1.23 for the = 1 and the flo = 0.2 
cosmological models, respectively, and we fix the shape pa- 
rameter of the CDM-like power spectrum (e.g., Bardeen 
et al. 1986) to F = 0.2. Figure 9 shows the compari- 
son between the predictions of PS and the observational 
mass functions. Observational data can be described by 
the PS models, except for the range of low-mass groups in 
the case oi = 1. Although it is well documented that 
the PS model which fits rich cluster data overpredicts the 
number density of low-mass halos compared to the sim- 
ulations (e.g.. Gross et al. 1998; Governato et al. 1999), 
the predicted difference is much smaller than the differ- 
ence between our observational MF and the (convolved) 
PS. Unless the observed difference is due to some problems 
of data incompleteness in the low-mass range, the case for 
the open model seems to be preferable. We emphasize that 
the comparison with PS predictions is not done here for 
determining the best fitting cosmological model. Instead 
it is aimed at verifying that our extension of the MF to 
group scales is reasonable. 

7 DISCUSSION 

7.1 Group Mass Determination 

We find that G93 groups can be described as systems 
in a phase of collapse, in agreement with previous results 
from earlier group samples identified in redshift galaxy cat- 
alogs (e.g., Giuricin et al. 1988; P92; Mamon 1994; Diafe- 
rio et al. 1993). We find that the presumably virialized 
region in groups is only <^ half of the radius sampled by 
galaxy data. 

In order to determine the dynamical status of groups 
and the corresponding mass correction we use the method 
of Giuricin et al. (1988), which has been also applied by 
P92. This method is based on simple numerical simula- 
tions, in which groups are not framed within a cosmologi- 
cal environment and the galaxies are represented by mass 
points starting from a spherical distribution with zero ve- 
locity (Giuricin et al. 1984). The comparison with ob- 
servations could be not straightforward since the observed 
galaxies may be affected by several environmental effects, 
e.g. tidal stripping, dynamical friction, and merging events 
which are not taken into account by these simple simu- 
lations. This is connected to the validity of the crucial 
hypothesis "mass distribution follows galaxy distribution" 
which is needed also for the standard virial mass estimate 
(see, e.g., Merritt 1987, 1988 for clusters) and which is also 
used by us in considering group masses contained within 
virialized regions (thus, our method of determining the dy- 
namical status is fully consistent with the mass estimate 
itself). 

Several N-body simulations of hierarchical clustering 
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have been performed in order to study galaxy systems in 
realistic situations. However, until recently, the poor res- 
olution of halos within dense environments has led to soft, 
diffuse halos that are rapidly dissolved by tidal forces (the 
so-called overmerging problem. White et al. 1987). This 
makes it difficult to compare simulations with observed 
group galaxies which are identified with halos. Recent 
simulations show that environmental effects can be impor- 
tant and affect the structure of individual galaxies (e.g., 
Moore et al. 1996; Tormen, Diaferio, & Sycr 1998; Colpi, 
Mayer, & Governato 1998). As for the global distribution 
of galaxies within clusters, Ghigna et al. (1998, 2000) per- 
formed very high-resolution simulations of clusters where 
the overmerging seems to be globally unimportant (the 
same seems to hold for small systems, see Moore et al. 
1999). Ghigna et al. foimd that the velocity dispersions 
of the halos agree with that of the dark matter particles. 
As for spatial biasing, they found that, at an early epoch 
of cluster formation, halos and dark matter have mimbcr 
density profiles of similar shape, while at the final time 
halos are anti-biased. In any case, since up to now these 
results concern virialized systems, studies of larger, unviri- 
alized regions should be awaited before reaching definitive 
conclusions for just forming groups. 

From the observational point of view, the luminosity 
segregation of galaxies in clusters, attributed to dynam- 
ical friction and merging, was found to concern only the 
brightest or very bright galaxies (Biviano et al. 1992; Stein 
et al. 1997). The overall global properties of clusters do 
not appear to depend on galaxy luminosities and also for 
groups virial masses are largely insensitive to luminosity 
weighting procedures (Giuricin et al. 1982; P92). 

On the other hand, there are several evidences that dif- 
ferent galaxy populations (e.g., early- and late-type galax- 
ies, blue and red galaxies, emission- or non-emission-line 
galaxies-hereafter ELGs and NELGs) show different spa- 
tial and velocity distributions (Biviano et al. 1997; Carl- 
berg et al. 1996, 1997b; den Hartog & Katgert 1996; 
Adami, Biviano, & Mazure 1998; Dressier et al. 1999; 
Girardi et al. 1996, 1998b; Mohr et al. 1996; Stein 1997; 
Koranyi & Geller 2000). As shown by the analyses of 
velocity dispersion profiles and spatial distribution, the 
galaxy component whose behavior most differs from the 
norm is that of very late-type galaxies or alternatively 
ELGs, being often ELGs very late types and vice-versa 
(Biviano et al. 1997; Adami et al. 1998). Biviano et al. 
(1997) suggested that the dynamical state of the ELGs re- 
flects the phase of galaxy infall rather than the virialized 
condition in the relaxed cluster core (cf. also Mahdavi et 
al. 1999). On the other hand, Carlberg et al. (1997a) sug- 
gested that, although differing in their distributions, both 
blue and red galaxies are in dynamical equilibrium with 
clusters, cf. also Mazure et al. (2000) who explain ELG 
dynamics by resorting to more radial orbits with respect 
to NELGs. 

As for global properties, as estimated within large sam- 
ples, using ENACS (ESO Nearby AbeU Clusters Survey) 
data, Biviano et al. (1997) found that the velocity dis- 



persion and virial masses based on ELGs arc, on aver- 
age, 20% and 50%, respectively, larger than those based 
on NELGs. However, due to the small fraction of ELGs 
(^ 10% of cluster members) the presence of ELGs does 
not strongly affects the estimate of velocity dispersion and 
mass of clusters. As for groups, we consider the 20 well- 
sampled groups analyzed by Mahdavi et al. (1999). They 
detected no velocity dispersion segregation between ELGs 
and NELGs together with a (not significant) decrease of 
masses by 20%, if ELGs are excluded from the sample. 
Therefore, even if global dynamical properties based on 
ELGs and NELGs can be significantly different, we ex- 
pect that the possible presence of ELGs in the groups we 
analyze hardly affects significantly, on average, our mass 
estimates. 

Finally, wc note that the dynamical status of the groups 
analyzed as well as the small number of galaxy members 
prevents us from applying refined analyses (used for clus- 
ters and well sampled groups), e.g. the determination of 
velocity anisotropics from velocity dispersion profiles and 
the Jeans equation. In fact, since the Jeans equation rigor- 
ously holds only in rc^gions being in dynamical equilibrium, 
this analysis was generally applied to galaxy clusters (e.g., 
Carlberg et al. 1996, 1997b; Ghardi et al. 1998b). This 
kind of analysis was also applied out to external regions 
of galaxy groups as a successful approximation (Mahdavi 
et al. 1999), but in any case it requires a large number of 
data. In this case Mahdavi et al. combined together data 
of a well-behaved subset of groups as selected by an analy- 
sis of the velocity dispersion profiles of individual groups. 
In this sense, the approach we use should be viewed as 
an alternative method for deriving masses in the case in 
which available data do not allow sophisticated analyses. 

1.2 Comparison with Previous MF 

At present, P92 is the only study of the mass dis- 
tribution function of loose groups. P92 determined the 
group mass function analyzing nearby groups {cz < 2000 
km s~^) of two catalogs based on the percolation method 
(38 and 21 groups by Geller & Huchra 1983 and Maia, 
da Costa, & Latham 1989, respectively) and one catalog 
based on the hierarchical method (107 groups by TuUy 
1987). The three catalogs have parent galaxy samples 
which considerably differ both for the selection criteria and 
the sky region covered. 

P92 noted significant differences between the group MFs 
resulting from the three catalogs of groups and claimed 
that these differences were not due to the inhomogene- 
ity of the catalogs, but rather to the choice of the group- 
selection algorithm, the percolation methods giving larger 
masses than the hierarchical one. 

Here, we consider the three group catalogs of G93, each 
having about 450-500 groups with cz < 5500 km s~^. For 
the determination of the MF, we retain only groups having 
\b\ > 20°, in order to avoid regions of high galactic extinc- 
tion. The homogeneity of the three catalogs, which come 
from the same parent galaxy sample, allows us to cast light 
on the compatibility of different selection algorithms. If we 
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take only the nearby groups [cz < 2000 km s~^, i.e. 117, 
124, and 136 groups) for a better comparison with P92, we 
find no diflference in MF among the three catalogs. This 
probably occurs because G93 chose the free parameters of 
the selection algorithms so as to obtain similar catalogs of 
groups. 

However, if we consider the whole catalogs (for which 
statistics is better, in particular at high masses) we find 
that percolation and hierarchical algorithms give really dif- 
ferent MFs, the former providing larger masses. A similar 
difference is also reported by P92 and it is proved here to 
be clearly due to differences in the algorithms. Indeed, it 
has been suggested that the drawback of percolation meth- 
ods is the inclusion in the catalogs of possible non-physical 
systems, like a long galaxy filament aligned close to the 
linc-of-sight, which give large mass estimates, while the 
drawback of hierarchical methods is the splitting of galaxy 
clusters into various subunits, which give small mass esti- 
mates (e.g., Gourgoulhon ct al. 1992). 

The difference among the three MFs is particular rel- 
evant in the high-mass range and leads to a fiattcr MF 
in the case of P groups. This effect is not seen in the 
analysis of P92 since they have no groups with M > 
4 X IQi^'/i^i Mq and less than 10 groups for M > 
lO^^h-^ Mq. 

As for the completeness, the inspection of the three 
differential MFs indicates that our samples are complete 
for M > 9 X 10^2^-1 Mq. Similarly, P92 assumed that 
the selection functions of both algorithms are efficient for 
1.1 X lO^-'/i-i Mq. 

As for the theoretical comparisons, P92 fitted their 
data to Press-Schechter predictions by assuming a unique 
power-law for the fluctuation power spectrum. However, 
since P92 did not give the amplitude of the fitted MF, a 
quantitative comparison with our results is not straight- 
forward. We only note here that their MF power-law slope 
lies in the range 1.5 ^ 1.7, being thus consistent with 
our result from the Schechter-like fit, 1.3^ a <; 1.6. As 
for the value of M*, P92 determined its value with quite 
large uncertainties, due to the small number of high-mass 
systems in their sample. 

8 SUMMARY AND CONCLUSIONS 

We analyze the three catalogs of nearby loose groups 
by Garcia (1993). This author identified groups in a 
magnitude-limited redshift galaxy catalog, which covers 
about ^ 2/3 of sky within cz = 5500 km s~^, by using 
two methods, a percolation and a hierarchical method. 
She tuned the free parameters of the group-selection al- 
gorithms in order to obtain two similar catalogs of groups 
and proposed a third catalog of groups defined as a com- 
bination of the two. Each catalog contains ~ 450-500 
groups. 

In agreement with previous works on earlier catalogs 
we find that these groups can be described as collapsing 

systems. Their typical sampled size is considerably larger 
than their expected virialized region. For all groups we 



compute the virial mass and correct this mass by taking 
into account the young dynamical status of these groups. 
We estimate corrected group masses, M, for two reference 
cosmological models, a flat one with = 1 for the mat- 
ter density parameter and an open one with flo = 0.2. 
For each of the three catalogs we calculate the group mass 
function, MF. 

Our main results arc the following: 

1. The number density of groups is not very sensitive 
to the choice of the group-identification algorithm. 
In fact, we find that the density of groups with 
M > 9 X 10^2 h^^ Mq, which is the adopted limit 
of sample completeness, ranges in the interval 1.3- 
1.9 X lO-^{h-'^Mpc)-^ for = I and it is about a 
factor of 15% lower for ~ 0.2. 

2. As for the MF shape, the percolation catalog gives 
a flatter MF than other catalogs. The difference 
decreases if we do not consider the most massive 
groups, for which reliable results come from galaxy 
cluster analysis. 

3. After obtaining the masses contained within the cen- 
tral, presumably virialized, group region by adopting 
a reduction in mass of ~ 30-40%, we do a com- 
parison with results coming from the virial anal- 
ysis of nearby rich galaxy clusters (Girardi et al. 
1998a). All three group MFs can be regarded as 
a smooth extrapolation of the cluster MF at M < 
4 X 10^'^h~^ Mq, which is the completeness limit of 
the cluster sample. Following a maximum-likelihood 
approach and taking into account mass uncertain- 
ties, we fit the Schechter expression for our MF on 
the whole mass range (groups with 9 x lO^^/i"^ Mq < 
M < 4x lO^"* h-'^ Mq and clusters with M > 4x 10^* 
h^^ Mq) and we obtain a slope of ~ —1.5 and a 
characteristic mass of ~ 3 x 10^** Mq. 

Our result strengthens the growing evidence in favor of 
the continuity of clustering properties from poor groups 
to very rich clusters (e.g., Ramella et al. 1989; Burns et 
al. 1996; Mulchaey & Zabludoff 1998; Ramella et al. 1999; 
Girardi ct al. 2000). 

Our resulting mass function of galaxy systems, which 
extends over two orders of magnitude, seems to be reason- 
ably described by Press-Schechter predictions of models 
which at larger masses describe the rich cluster mass func- 
tion (in particular, in the case of the open model). 

The analysis of wide forthcoming group catalogs, e.g. 
the UZC (Updated Zwicky Catalog) groups by Pisani 
et al. (2000) and the NOG (Nearby Optical Galaxies) 
groups by Giuricin et al. (2000), as well as extensive spec- 
troscopy observations of individual groups (e.g., Zabludoff 
& Mulchaey 1998; Mahdavi et al. 1999), will be of great 
aid in improving the determination of the group MF, in 
particular at the low-mass range which is affected by larger 
uncertainties. 

Furthermore, future studies of N-body cosmological sim- 
ulations for different cosmological scenarios will improve 
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the understanding of non-virialization effects on the esti- 
mate of group masses (e.g., Diaferio et al. 1999). 
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